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We study controllable exchange coupling between two singlet-triplet qubits. We start from the 
original second quantized Hamiltonian of a quadruple quantum dot system, and obtain the effective 
spin-spin interaction between the two qubits using the projection operator method. Under a strong 
uniform external magnetic field and an inhomogeneous local micro-magnetic field, the effective 
interqubit coupling is of the Ising type, and the coupling strength can be expressed in terms of 
quantum dot parameters. Finally, we discuss how to generate various two-qubit operations using 
this controllable coupling, such as entanglement generation, and a controlled-NOT gate. 

PACS numbers: 03.67.Lx, 73.21. La, 75.10.Jm 



I. INTRODUCTION 

A quantum computer is more efficient than a classi- 
cal computer in solving certain problems such as prime 
factorization^ A variety of physical systems have been 
suggested as qubits, the building block of a quan- 
tum computer^ Spin qubits have attracted wide atten- 
tion for more than a decade 3 - - — because of the well- 
developed spin resonance techniques^ for coherent control 
and the strong exchange coupling between electron spins 
in semiconductors £ In recent years, there have been a 
multitude of experiments verifying spin control^— spin 
coherence ^ 1 1 ' 14 ' 17 ' 18 and spin measurement i 19 i 20 

Singlet-triplet qubit (STq qubit) is a logical qubit en- 
coded in the two-spin singlet ^j(|t4-) — lit)) an d unpo- 
larized triplet -^jdtl) + lit)) states in a double quan- 
tum dot (DQD)i 21 ' 22 Since both states have zero mag- 
netic quantum number, this qubit is insensitive to noises 
in magnetic fields that are uniform or slowly varying in 
space. More importantly, this qubit can be initiated 
with high fidelity and manipulated with prccisioni 10 i 13 
Its fast measurement has been demonstrated as wel h 10 i 19 
Coupling of two singlet-triplet qubits can be imple- 
mented capacitively^ 3 - or through exchanged Recent 
theoretical 2 ^— and experimental^ explorations are fo- 
cused on the capacitive coupling. However, this approach 
has a sensitive inter-dependence of the two-qubit cou- 
pling strength and the single-qubit dephasing rate ) 16 i 26 
with stronger coupling leading to faster dephasing and 
vice versa. To overcome this problem, additional as- 
sistances such as a long-distance coupler have been ex- 
plored to provide enhanced capacitive coupling without 
increased dephasing 

The difficulty in achieving strong capacitive coupling 
without increased decoherence prompts us to revisit the 
coupling scheme in the original proposal, using the ex- 
change interaction between the two STq qubits. Here we 
address one of the main issues raised in the original pro- 
posal about the exchange coupling, i.e., the issue that 
it potentially leads to leakage out of the STq qubit sub- 



space. Specifically, we calculate the effective coupling 
between two STq qubits, and identify the conditions un- 
der which the qubit leakage can be minimized. We also 
construct several useful two-qubit gates based on the ex- 
change interaction. 

The paper is organized as follows. In Sec.|Hl we derive 
the exchange coupling between two electron spins from 
the second quantized Hamiltonian of a DQD system. In 
Sec. IIII1 we extract a controllable coupling between two 
STq qubits from an exchange-coupled quadruple quan- 
tum dot system. In Sec. IIV1 we study how to generate 
entanglement and build controlled-NOT gate using this 
controllable coupling. At last, we give a summary in 
Sec. El 



II. SINGLET-TRIPLET QUBIT IN A DOUBLE 
QUANTUM DOT 

The two-electron spin Hamiltonian has been studied 
extensively in the past decade^i 2 ^— and its form in the 
STq subspace has also been explored^ 2 - Here we quickly 
summarize the results as a starting point for our calcula- 
tions in the next section and derive an analytical expres- 
sion for the two-electron exchange coupling that includes 
all the tunable parameters. 

As discussed in Appendix [AJ we derive the two- 
electron effective Hamiltonian from a generalized Hub- 
bard model<22r— In the so-called (11) regime, where each 
quantum dot contains one and only one electron, the ef- 
fective spin interaction Hamiltonian for two electrons in 
a DQD takes the simple Hciscnberg exchange form, 

HcS = JSi ■ S2, (1) 

with the exchange coupling between the two electrons 
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Here t is the single-electron tunneling across the DQD, 
J t is the single-electron tunneling in the presence of a 
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second electron, U is the on-site (i.e., intradot) Coulomb 
repulsion, U' is the interdot Coulomb repulsion, Ae is the 
single electron ground-orbital energy difference between 
the two dots (interdot bias) , and J e is the direct exchange 
interaction of the two electrons across the DQD. 

As discussed in Appendix [3] the exchange splitting 
J in Eq. @ is derived from the generalized Hubbard 
model and is quite complete; it includes the effects of 
interdot tunneling, interdot bias, and both on-site and 
off-site Coulomb interactions. The expression for J con- 
tains two parts. The first part, 4(i- J t ) 2 /(U- U' - \ Ae\), 
is the antiferromagnetic superexchange between the two 
dots. The second part, — 2J e , is the ferromagnetic direct 
exchange between the two electrons from their Coulomb 
interactions. The value of total exchange J can be ei- 
ther positive or negative, depending on the values of the 
parameters t, J t , Ae, J e , U, and U' . This switch in the 
exchange coupling between anti-ferromagnetic and fer- 
romagnetic interactions has been previously observed in 
calculations based on various levels of molecular orbital 



approximations 
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It is important to note here that 



the zero point for the interdot bias here differs from the 
convention in the singlet-triplet qubit community. Here 
Ae = means that the single-electron ground states are 
on resonance across the two dots, while in the conven- 
tional definition the zero point is where the (11) singlet 
and the (02) singlet states have the same energy when 
tunneling is neglected: el + er + U' = 2er + U. Thus 
the two bias definitions arc shifted by U — U'. 

Including the external magnetic field (applied along 
the z direction) and the Overhauser field from the lattice 
nuclear spins, the complete two-spin Hamiltonian is 

H = j e BSl + le BS z 2 + JSi ■ S 2 + Si ■ hi + S 2 ■ h 2 , (3) 

where 7 e is the electron gyromagnetic ratio, and h 1( - 2 ) 



/,l(2) 1 /,l(2) 



is the nuclear field operator in each dot, 



with Ai t i(2) being the hyperfinc coupling between each 
electron spin and its nuclear spin bath. 

When the applied magnetic field is large, ^ e B S> 
J,|hi(2)|, the singlet S and the unpolarized triplet Tq 
states are isolated from the polarized triplet states T±. 
Under this condition, S and To can be used as the two 
basis states to encode a single logical qubit, i.e., the STq 
qubit. Within the STq subspacc, the qubit dynamics is 
governed by the Hamiltonian 3 - 2 - 



Sh z 



(4) 



where r z = \Tq)(T \-\S)(S\ and r x = |S)(T | + |T )(S| 
are the Pauli operators of the STo qubit, with X 1 Y and 
Z representing the three axes of the STq Bloch sphere. 
6 h z = (h\ — is the Overhauser field difference (in 

general it is the total magnetic field difference that could 
come from both the Overhauser field and an externally 
applied inhomogeneous field) across the DQD along the £ 
direction in real space. The exchange splitting J and the 
magnetic field difference 8 h z provide universal control 
over a single STq qubit. 




FIG. 1: (Color online) Four quantum dots are designed to be 
coupled as an array, where only single orbital level is consid- 
ered in each quantum dot, and on-site (i.e., intradot) Coulomb 
repulsion is denoted by U. Here we allow the tunneling cou- 
pling £23 between dots 2 and 3, which is controlled by the 
external gate voltage, i.e., the height of potential between 
dots 2 and 3. 



III. CONTROLLABLE COUPLING OF TWO 
SINGLET-TRIPLET QUBITS 



We now study the exchange coupling of two STq 
qubits. The system we consider is a linearly coupled 
quadruple quantum dot (see Fig. [TJ, with dots 1 and 
2 encoding the first qubit, and dots 3 and 4 encoding the 
second. In addition to the tunneling between dots 1 and 
2, and the tunneling between dots 3 and 4, we further al- 
low tunnel coupling between dots 2 and 3, so the two STq 
qubits are now coupled. This is the main difference be- 
tween this coupling method and the previous capacitive 
coupling scheme! 16 ' 23 " — 

As discussed in Appendix [51 the effective interaction 
Hamiltonian for the quadruple quantum dot system is 
that for a linear Heisenberg spin chain of four nodes: 



J12S1 • S 2 + J 2 3S 2 ■ S3 + J34S3 • s 
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where the exchange couplings are 
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Here U is the intradot Coulomb repulsion. The in- 
terdot Coulomb repulsions are U12 = U34 = U', and 
U23 = U". The parameters t k ,k+i, J? M1 , J^ k+1 , 
and Aefc,fc+i characterize the nearest-neighbour tunnel- 
coupling, occupation-modulated tunneling, direct ex- 
change, and interdot bias, respectively. The exchange 
couplings Jk,k+i are tunable via tk,k+i and Ae^^+i, 
while the Coulomb parameters U, U', U", J e , and J t 
cannot be tuned easily (see Figs. [T] and E]). 

After including all the magnetic interaction terms, the 
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total four-spin Hamiltonian reads 

2 4 

k=l fc=3 
+ J12S1 • S2 + J23S2 ' S3 + J34S3 • S4 

4 

+ ^S,-h fc , (7) 

fe=i 

where hfc is the nuclear Overhauscr field in dot k. We 
have also introduced a local magnetic field for each STo 
qubit, so spins 1 and 2 are in a local micromagnetic field 
B m , while spins 3 and 4 are in an opposite micromag- 
netic field — B m , in addition to the overall uniform mag- 
netic field B. This local field can be produced by local 
micromagnetS ) 12 i 36 ' 37 and prevents qubit leakage, as dis- 
cussed below. 

In general, the Hilbcrt space for the four electron spins 
has a dimension of 16. When a strong uniform external 
magnetic field B along the z direction is applied, i.e., 
-f e B S> |hfe| and j e B Jk,k+i, spaces with different 
total magnetic quantum number S z = X/fc=i ^fc are de- 
coupled. The ST qubits satisfy S z — 0, so we focus 
here on the S z = sub-Hilbert-spacc. It is spanned by 
\SS), \ST ), \T S), \T T ), |T + T_>, and |T_T+). No- 
tice that the two-ST -qubit Hilbert space, spanned by 
\SS), \ST ), \T S), \T T ), has only four dimensions and 
is only a subspace of the S z — space. Hamiltonian 
([5]) shows that states in the STo subspace are coupled to 
|T+T_) and |T_T+) by J 23 S 2 • S 3 . To prevent leakage 
into these two states, we introduce the local magnetic 
field ±B m to separate them energetically from the two- 
qubit Hilbcrt space. When J23 <C 2j e B m < 2j e B, the 
ST states would be decoupled from |T + T_) and |T_T+), 
so two STo qubits would evolve without loss when they 
are coupled through the ^-coupling. 

Projecting Hamiltonain ([7]) into the two-qubit Hilbert 
space, we obtain the following effective Hamiltonian 
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« le {B + B m ) S z k + le{B - B m ) S z k 

fc=l k=3 
+J12S1 ■ S2 + T34S3 • S4 + J23S2SI 
4 

+ ^S fe -h fc , (8) 
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where J^ 3 



J23 + (J 2 3) 2 /( 4 7eS m ) w J 23 - The qubits arc 
now well defined even when they are coupled, with spins 
1 and 2 forming one STo qubit (denoted as a), spins 3 
and 4 forming another STo qubit (denoted as 6), and the 
interaction T^Sf Sf couples these two STo qubits. Fig- 
ure [5] shows a diagrammatic representation of this pro- 
jection process. The four linearly-coupled electron spins 
in the specially designed magnetic field [Fig. [2ja)] are 
equivalent to two coupled pseudo spins [Fig. 0[b)], with 
each pseudo spin representing an STo qubit. Expressing 
all the electron spin operators of Eq. ([5]) in terms of the 



S 1 s 2 

(b) 



- b 4 
^ f-^-^ 



B 



()■ 



J23 X aTb/4 



.3 04 



■o 



FIG. 2: (Color online) (a) Exchange-coupled four electron 
spins derived from a linearly coupled quadruple quantum dot 
system. The exchange couplings Jk,k+i are controllable by 
external gate voltages. A strong uniform external magnetic 
field B is applied along the z direction. Also, we apply lo- 
cal micro-magnetic field to each spin; spins 1 and 2 are in a 
local micro- magnetic field B m , while spins 3 and 4 are in an 
opposite micro-magnetic field — B m . (b) Under this carefully 
designed magnetic field, effectively coupled two pseudo spins 
(each pseudo spin represents a STo qubit) are extracted from 
the linearly coupled four spins shown in (a). The pseudo spin 
T a is formed by Si and S 2 , and 1% by S3 and S4. The cou- 
pling between these two pseudo spins is generated from the 
exchange coupling J23 between spins 2 and 3, and the resulted 
controllable coupling is of the Ising XX type. 



pseudo-spin operators, we obtain the Hamiltonian of two 
coupled STq qubits 
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This is the Hamiltonian for a transverse field Ising model. 
In other words, tunnel-coupling of two STo qubits is 
equivalent to a controllable effective quantum Ising cou- 
pling between them. 

For a two-spin qubit, the singlet and triplet states are 
not the only way to encode a logical qubit. An alterna- 
tive is to encode the logical qubit in the |tl) and |4t) 
states. This encoding is less susceptible to charge noise, 
since the two-electron orbital wave functions and the as- 
sociated charge distributions are identical for 1 7 J,) and 
states, but differs for S and To states. Indeed, the 
recent experiment that demonstrated a 200-/xs coherence 
time in a GaAs double dot is for these two states*^ For 
this encoding scheme, the effective coupling Hamiltonian 
takes the form 



tt _ J 12 X , J 3i X 
tl — ——T„ +— — T h - 



J. 
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+6h z b r b z . (10) 



It is still the transverse field Ising model Hamiltonian. 
Thus the pulse sequences for two-qubit gates would be 
similar for this encoding as the STo encoding. 

Notice that the coupling we have here has a form that 
differs from the capacitively coupled STq qubits, where 
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the total Hamiltonian corresponds to a longitudinal field 
Ising model (or classical Ising coupling)! 16 ' 26 

In our exchange coupling scheme, the coupling strength 
J23/4 is only limited by the external fields J23 <C 
2~f e B m < 2^/ e B, but not by the single-qubit level spacings 
J12 and J34. Also, J23 can be tuned to be comparable to, 
or even larger than, J12 and J34. Indeed, J12 and J34 are 
allowed to vanish altogether J 12 = J34 = 0, and the sys- 
tem will still remain in the STq qubit space. This is in 
strong contrast to the capacitive coupling scheme^ 3 . - — 
where the coupling strength is limited by J 12 and J34, 
with J23/J12, J23/J34 ~ 10~ 2 in a recent experimental 
demonstration (see Ref . Htl ) . 

To estimate the maximum value of the achievable 
coupling strength of the interqubit coupling ^-t^t^ , 
we take a GaAs quantum dot structure as an exam- 
ple. With current experimental technology, a magnetic 
field gradient of ~ 20 — 30 mT can be realized . 12 i 36 ^I 
The corresponding Zeeman energy gradient is 2j e B m « 
120 — 180 MHz, so the maximum coupling strength is 
about J23/4 w 10 MHz. In principle, the magnetic field 
gradient 2"/ e B m can be several times larger, so an in- 
terqubit coupling close to 100 MHz (0.4 fieV) should be 
possible. 

IV. GENERATING TWO-QUBIT GATES 

Two-qubit gates are essential for universal quantum 
computing. In this section, we discuss how to use the 
Ising interaction derived in the previous section to gen- 
erate Bell states and to construct a controllcd-NOT gate. 

A. Entanglement generation via free evolution 



otherwise. Such a process can be realized in a quantum 
dot device by tuning the tunneling coupling i 2 3 between 
dots 2 and 3 in Eq. ((6]). If the two qubits are initialized 
to a desired product state |*(0)) = |ta-U)r^^ after the 
r-period evolution, the final state is 

\*(T))=e- iHT \t a U) = cos^\taU)+ism^\UW- 

(13) 

This is generally an entangled state for the two STq 
qubits. In particular, when J'r cn t/4 = nir ± 7r/4, two 
maximally-entangled Bell states are obtained: 

^(|T«4*> ± i\iaW) = e- iHT "*\tah). (14) 

If the initial state is set to |\P(0)) = |tat&)i f ree evo- 
lution under Hamiltonian (jlip leads to entangled states 
between | tatfc) an d l-J-oJ-b)- Specifically, when £aT en t = 
mr + 7r/2, we obtain the following entangled state: 

cos 0\ tatb) + sin 0\ UU) = e~ iHT ^ I t.tfc) • (15) 

The degree of entanglement in this case depends on the 
ratio J' j J. Since J' can be even larger than J, we can 
imagine tuning of the couplings to J' = 4J, such that 
6 = — 7r/4. Now the final state is another maximally- 
entangled Bell state. 

In short, starting from different unentanglcd initial 
states, a pulsed ./'-gate can generate different maximally- 
entangled Bell states. Since J' is limited by the interdot 
magnetic field gradient, it may not be feasible to reach 
the J' > J regime in GaAs while neglecting the local 
nuclear fields 8h z . However, this limit should be much 
easier to achieve in Si, where the nuclear Overhauser field 
is much smalleri 14 i 39 



For entanglement generation, we tune the qubit split- 
tings to J 12 = J34 = J 3> 5h z (putting the two qubits in 
resonance after neglecting the nuclear fields) and write 
J23 = J' ■ The two-qubit Hamiltonian now takes the sim- 
pler form 

11=^ + 1^-^^. (11) 

The eigenfunctions and the corresponding eigenvalues of 
this Hamiltonian are 



H) = cos-|t„t6>+Bk-|UW, Ei = VJ 2 + U74) 2 , 
1 2) = sin °- 1 fat*) - cos °- 1 UU) , E 2 = - VJ 2 + ( J'/4) 2 , 
|3> = ;^(|UU> - U«tb», E3 = \f, 
|4> = ^(|UU> + Uatb>), E i = ~J', (12) 

where 9 = arctan(— J'/4J). Since J' is controllable, we 
can turn it on for a period of time r, and keep it off 



B. Controlled-Not gate 

One of the most commonly used building block for uni- 
versal quantum circuits is the controlled-NOT (CNOT) 
gate^°- which can be generated relatively straightfor- 
wardly using an Ising-type interaction^ In our case, 
we can adopt two different approaches to construct a 
CNOT gate. The first approach starts with Eq. (jTO . 
After two single-qubit rotations U\ — e 4 "f( T a+ T b) anc l 
U 2 = e -^( T ° +T b\ Hamiltonian dTTJ) can be transformed 
to H' 

e -iH'r = U^e-iHTjj^ (lg) 

with 

H t =-±T*-ii?-£tt. (17) 

The pulse sequence of a CNOT gate using this Hamilto- 
nian is well know n 41 ' 42 

^CNOT = Hbe- 1 ^ e 1 ^ e- tH ' T e- l7TT ^ 2 e- lH ' T H b , (18) 
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where H}, = (r b + ) is the Hadamard gate on qubit 
b, the duration of the two-qubit gate is r = 7r/(4J) x 
*/ (4n) 2 — (2m — l) 2 , where m and n are integers, and 
(j) = (7r/4)(2m — 1). The coupling parameters must be 
tuned to the ratio of J/ J' = y/[4n/(2m- l)] 2 - 1/4. No- 
tice that the two-qubit gate is employed twice in this 
pulse sequence. 

Alternatively, we can tune the parameter J in Eq. 
to J = 0, which makes the construction of a CNOT gate 
rather simpler. The local nuclear Overhauser field 5 h z 
can be completely eliminated by a Hahn spin echo4^ The 
spin echo (the pulses are simultaneously applied to the 
two ST qubits respectively) under Hamiltonian ^ leads 
to a two-qubit quantum gate e*T T a T b ; where we have 
chosen J'r/4 = 7r/4. The CNOT gate can then be real- 
ized as follows 



jjab 



CNOT 



(19) 



Note that the two-qubit gate is applied only once, in 
contrast to our first approach (using a transverse field 
Ising Hamiltonian), or CNOT gates generated by using 
other types of exchange interactioni 8 i 44 i 45 

Other important two-qubit gates can be achieved by 
combining the controllcd-NOT gate and single-qubit 
gates. For example, the controlled- Z gate is given 
by = H b U^ NOT H b ^^. and the swap gate is 

a combination of three controllcd-NOT gate: US^ = 

jjab Tjba J job 40 
17 CNOT ^CNOT u CNOT - — 



DISCUSSIONS AND SUMMARY 



more complete formula for the exchange coupling be- 
tween two confined electrons in a DQD. The tunable 
exchange coupling depends explicitly on intradot and in- 
terdot Coulomb repulsions, tunnel couplings, and inter- 
dot bias. These two exchange-coupled electron spins can 
then be used to encode a singlet-triplet qubit in a strong 
magnetic field. We then use the generalized Hubbard 
model to construct the spin Hamiltonian for two tunnel- 
coupled STq qubits, and derive a general expression for 
the exchange couplings between neighboring dots in a lin- 
ear four-dot configuration. This exchange-based coupling 
scheme leads to an effective Ising interaction between the 
two STq qubits. The coupling strength is limited by the 
magnetic field gradient between the two DQDs, but not 
limited by the single-qubit energy splittings, and we es- 
timate that the interqubit coupling strength can reach 
the range of 100 MHz. We also explore how to generate 
various two-qubit gates with this coupling. In particular, 
we discuss how to construct the controllcd-NOT gate and 
how to generate Bell states. 

The results we obtain in our study of a four-quantum- 
dot system indicate that exchange coupling is a viable 
alternative to capacitive coupling for STq qubit. It is 
more tunable, and its coupling strength is not directly 
connected to single-qubit decoherence. With these two 
favorable features, this exchange-based coupling scheme 
for the STq qubits deserves closer experimental scrutiny. 

Note added: As we were finalizing our manuscript, we 
became aware of the interesting work by Klinovaja et 
al.r^ which also dealt with the exchange coupling of STq 
qubits, and included the effect of spin-orbit coupling as 
well. 



In our study, we did not consider qubit decoherence 
due to the environment. In real quantum dot devices, the 
main decoherence sources are the lattice nuclear spins 4 ^ 
and the background charge fluctuations 4 ' 32 ' 48 For a DQD 
with a finite exchange coupling </, charge noise is partic- 
ularly harmful. In our coupling scheme, however, we can 
tune the system to the small J limit in order to minimize 
this decoherence effect. Furthermore, it may be possible 
to remove/suppress decoherence with various dynamical 
decoupling schemes applied to each qubit (see, e.g., Refs. 
I20I and [49T ). In short, it is reasonable to expect that 
this exchange-based qubit coupling scheme should lead 
to better decoherence performance for the STq qubits 
than the capacitive coupling scheme. 

The transverse field Ising model has been exten- 
sively explored in the areas of adiabatic quantum 
computing^ - — and quantum phase transition.— Our re- 
sults here show that the transverse field Ising model can 
be realized in a tunnel-coupled STq qubit system, making 
it a potential platform or quantum simulator for imple- 
menting adiabatic quantum computing and investigating 
quantum phase transitions. 

In summary, starting from the recently introduced gen- 
eralized Hubbard model^^r— we analytically derive a 
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Appendix A: Derivation of the effective interaction 
Hamiltonian for two electron spins 

In this Appendix we derive the effective spin Hamilto- 
nian for the low-energy dynamics of two electrons con- 
fined in a semiconductor DQD potential V(r), as shown 
in Fig. [31 The two electrons interact with each other via 
the Coulomb interaction and can be used to achieve the 
so-called STq qubit. The effective mass Hamiltonian for 
such a two-electron system takes the form 

H = E A + E + 1 fi2 , . ( A1 ) 

^ 2m* ^ e s ri - r 2 

1=1 e 1=1 s| 1 z ' 
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bard model^^i 



FIG. 3: (Color online) Double quantum dot model. Two 
electrons are confined by a two- well potential V(r), where 
each well represents a quantum dot and a*, is the center of 
the dot k. Here only a single electron ground-orbital energy 
level is considered for each dot, and t represents the tunneling 
coupling between the two dots, which is controllable by tuning 
the potential height between the two dots. U is the intradot 
Coulomb interaction, and As is the single electron ground- 
orbital energy difference between the two dots. 



where m* is the effective mass of the conduction electrons 
and e s is the static dielectric constant of the material. 
Since the total spin is a good quantum number for this 
Hamiltonian, the two-electron states can be classified into 
singlets and triplets, with symmetric and anti-symmetric 
orbital wave functions. The problem can then be solved 
numerically with the molecular orbit method (or ana- 
lytically within the Heitler-London and Hund-Mullikan 
approximations) .— ~— The splitting between the ground 
singlet and triplet states is the exchange coupling be- 
tween the two confined electrons. With large on-site 
Coulomb interaction and single-particle excitation en- 
ergy, the problem can also be solved using the Hubbard 
model, starting with a second quantized Hamiltonian by 
defining a field operator <J(r) = Ylk=i,<? c fc^<?K r ~ a fe)x<r, 
where </>(r — a&) is the single electron ground orbital wave 
function in dot k centered at afc ) 29 i 56 and \a is t ne s P m 
cigenfunction that describes the electron spin degree of 
freedom. While in this model higher-energy orbital states 
are neglected, it does have the advantage that it can be 
easily extended to larger multiple-dot systems. The sec- 
ond quantized Hamiltonian reads 



H 



2m* 



+ V(r) tf (r) + 



- f drdr'^(r)^(r')—^ -*(r')*(r). (A2) 



Expressing the Hamiltonian H with the fermionic cre- 
ation and annihilation operators for each single-dot 
single-electron eigenstate, one obtains a generalized Hub- 



H = H e + H t + Hu + Hj, 

2 



H 



fc=l,, 



ff t = t5^(4 r c to +H.c.) J 

cr 

Hu = Unifmi + Un 2 ^n 2 i + U\nv\ + nii)(n 2 -\ + "24.), 
Hj = - J e (nifn 2 f + 7114.7124.) - (J e c}^4t c 24Cit 

+ J P 4 t 4_|_citci4 + ^2 Jtn ia c\ s C2s + H.c.) . (A3) 

icj 

Here H e and H t represent respectively the single electron 
energy level of each dot and the tunnel-coupling between 
the two dots. Together they make up the single-particle 
part of the Hamiltonian. The two-body part, i.e., the 
Coulomb interaction, is described by Hu and Hj. Pre- 
vious studies often started from H = H e + H t + Hu, 
which only gives an anti-ferromagnetic exchange coupling 
between the two electrons.— However, numerical results 
based on either Heitler-London approximation or Hund- 
Mulliken molecular-orbital method have shown that the 
exchange coupling can be ferromagnetic in certain pa- 
rameter regimes.— _ — Here, starting from the general sec- 



ond quantized Hamiltonian in Eq. (|A3[) . we can obtain 
a more precise expression for the exchange coupling be- 
tween the two electrons. 

In the limit of (U — U') ^> t, there is no charge degree 
of freedom for the two electrons; each dot confines one 
and only one electron (i.e., double-occupation in each dot 
is forbidden), and the two electrons are only allowed in 
the (11) charge configuration. Below we derive the effec- 
tive spin-spin interaction between the two electrons, i.e., 
the exchange coupling, by using a projection operator 
method. We define a projection operator 



P = JJ["fct(l - n n) + n ki {\ - Tifet)], 



(A4) 



fe=i 



which projects the Hilbert space of the two electrons to 
the subspace H. for (11) charge configuration. Using a 
standard procedure^ we derive an effective spin Hamil- 
tonian that describes the dynamics of TL, 

H eS = PHP - PHQ(QHQ - E)~ r QHP, (A5) 

where Q — 1 — P. After some algebra, we obtain the 
following relations: 



PHP = J2 e k- 

k=l 

PHQ = PH' t Q, 



U' - i J e - 2J e (S! • S 2 + i), 
QHP = QH' t P, 



(A6) 



(t ~ Jt) J2<t( c i<j c 2<? + H.c). In other words, 
under projection operations, the Jt term in Hj has the 



where H' t 
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same effect as H t . In addition, under the lowest order 
approximation, QHQ — E can be replaced by the differ- 
ence between the first charge excitation energy 2ei + U 
and the ground state energy e\ + £2 + U', where £2 > £1. 
Thus, we have 



QHQ-Ek U-U'- \Ae\, 



(A7) 



where Ae — £2 — £1. Neglecting the constant terms, we 
obtain the effective Hamiltonian 

„ _ 2 r S s (PH' t Q)(QH' t P) 

H cff - -2J e bi • S 2 - u _ u , _ | Ag | ■ ( A8 ) 

Using the properties Q 2 = Q, Q = 1 -P, and PH[P = 0, 
it follows that 



H e ff = — 2J e Si • S 2 — P 



U-U'- |Ae 



P. (A9) 



Now we focus on simplifying the second term in the above 
equation, 

PH' 2 P = (t- J t f [Pclc la/ PPc 2a cl,P 

<J(j' 

+Pc\ a C 2 a' PPc X ac\ a , P] . ( A10) 

We can use the following identities: 

Ck a c\ a , = 5aa> [l - X ("fet + n *i)] _ S fe ' > (All) 

where the spin operator S fc = | ^ CT (T / c^ <T <r - -/Cfc <7 /, with 
<t being the Pauli operator, and <t cto -/ its matrix elements. 
The number operators satisfy PQ2 a n k<x)P = 1- After 
some algebra, we obtain 



4(t- J t ) 2 (i-S 1 -S 2 ). 



(A12) 



Substituting this equation into Eq. (|A9|) . and neglecting 
the constant terms, we finally obtain the effective inter- 
action Hamiltonian for the two electron spins 



where 



Heff — </Si • S2, 



U-U'- |Ae| e 



(A13) 



(A14) 



is the exchange splitting between the two electrons. 
These are the results given in Eqs. ([TJ and 



Appendix B: Derivation of the effective interaction 
Hamiltonian for four electron spins 

As discussed in the main text, the coupled qubit sys- 
tem that we consider is a linearly coupled quadruple 
quantum dot (see Fig. [1]), with dots 1 and 2 encoding 
the first qubit, and dots 3 and 4 encoding the second. 
In addition to the coupling between dots 1 and 2, and 
that between 3 and 4, we further allow tunnel coupling 
between dots 2 and 3, so the two STq qubits are tunnel 
coupled. 

Similar to the DQD case discussed in Appendix |A1 we 
define a field operator *(r) = £)fc=i,o- c k>a <j>(r - a k )xa 
by using the ground orbital state <p(r — a k ) in each dot. 
The second quantized Hamiltonian of this quadruple dot 
system is then 

4 3 
H = £kc\ a c ka + ^ t k ,k+i(c{ a c k+h(T + H..c.) 

k—l,a k— 1,(7 

fe=l k=l 
3 

+ y^^fc.fc+i("fet + n ki)(n k +i,t + n k+i,i)- (Bl) 

k=l 

Here e k again describes the single electron ground-orbital 
energy in dot k, t k k+ i is the tunneling coupling between 
two neighboring dots k and k + 1, U is the intradot 
Coulomb repulsion, and both Hj' k+1 7 which is defined 
in Eq. (|A3[) . and U k ,k+i describe the Coulomb interac- 
tions between two nearest-neighbor dots k and fc + 1. For 
simplicity, we consider the symmetric case (between the 
two STq qubits) with the parameters U12 — U34 = V, 

TJ — TJ" T 12 — 734 t QT1 j 723 _ 77 

'-'23 — u 1 J e,p,t ~ J e,p,t ~ J e,p,ti <" lu J e,p,t ~ J e,p,t' 

Under the strong Coulomb interaction conditions, with 
(U-U') > ti2,t 34 and (U-U") > t 2 3, each dot confines 
only one electron J2a n k°- = 1- Other charge configura- 
tions are not allowed, and these four electrons have only 
the spin degrees of freedom. 

We now apply the same projecting procedure as for 
a DQD in Appendix \K\ in order to obtain the effective 
spin Hamiltonian for the four electrons. The projection 
operator is defined as 



P = [iKtt 1 _ n u) + n-ni 1 - n kt)}- 



(B2) 



k=l 



After applying this operator, the Hilbert space for the 
four electrons is reduced to one involving only the (1111) 
charge configuration. Following the same procedure as in 
Appendix El we obtain the effective interaction Hamil- 
tonian for the four electron spins 

-ffeff = JiaSi • S2 + J23S2 • S3 + J34S3 • S4, (B3) 
where 



Jk.k+l = 



4{t k>k+1 - J, 



fc,fc+l\2 



U -U' + U"- I Asu k+1 1 



2J; 



k,k+l 



(B4) 
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These are the results given in Eq. (J5)) and ([6]). The pa- 
rameters ifcfc+i, J t ' , Jg ,fe+1 , and Aefc fc +1 characterize 
the tunneling between nearest-neighbor dots, occupation- 
modulated tunneling, spin exchange, and single electron 
ground-orbital energy difference, respectively. We em- 
phasize here that all the exchange splittings Jk.k+i arc 



controllable. On the other hand, once the quantum dot 
device is designed, the Coulomb-interaction parameters 
U, U' , U", J e , and Jt are not easily tunable. The directly 
tunable parameters tk,k+i an d Aei^k+i are controlled by 
the external gate voltages (see Figs. [T]and[3]). 
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